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PEP-II asymmetric-energy collider. The data were recorded with the BABAR detector at center- 
of-mass energies 10.58 and 10.54 GeV, and correspond to an integrated luminosity of 454 fb -1 . 
We investigate the J/ip-K + ir~ mass distribution in the region from 3.5 to 5.5 GeV/c 2 . Below 3.7 

GeV/c 2 the tp(2S) signal dominates, and above 4 GeV/c 2 there is a significant peak due to the 
Y(4260). A fit to the data in the range 3.74 - 5.50 GeV/c 2 yields a mass value 4244 ± 5 (stat) 

± 4 (syst)MeV/c 2 and a width value 114 t\t (stat) ± 7(syst) MeV for this state. We do not 
confirm the report from the Belle collaboration of a broad structure at 4.01 GeV/c 2 . In addition, 
we investigate the ir + ir~ system which results from Y(4260) decay. 
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1 The observation of the X(3872) [TJ, followed by the dis- 21 

2 covery of other states such as the Xc2(2-P)(3930) [2], the 22 

3 V(3940) [3 , and the A" (3940) 0], has reopened interest 23 

4 in charmonium spectroscopy. These resonances cannot 24 

5 be fully explained by a simple charmonium model [5]. 25 
e The V(4260) was discovered [B] in the initial-state- 26 
7 radiation (ISR) process e+e~ -> 7i SR V(4260), V(4260)^ 27 
a J/ipTr + Tr~ . Since it is produced directly in e + e _ annihi- 28 

9 lation, it has J PC = 1 . The observation of the decay 29 

10 mode J/ipTT°TT° [7] established that it has zero isospin. 30 

11 However it is not observed to decay to D* D* [B , nor 31 

12 to D* s D* s 9 J , so that its properties do not lend them- 32 

13 selves to a simple charmonium interpretation, and its 33 

14 nature is still unclear. Other interpretations, such as 34 

15 a four-quark state [TD1 QTJ , a baryonium state [T2] , or 35 

16 a hybrid state [13], have been proposed. However if 36 

17 the Y(4260) is a four-quark state it is expected to de- 37 
is cay to DfDj [TT], but this has not been observed [9]. 38 

19 An analysis of the reaction e + e~ — > J/i/jtt + tt~ [14] 39 

20 which confirms the Y(4260), suggests the existence of 40 



a broad state with mass m = 4008 ± 40 MeV/c 2 
and width T = 226 ± 44 ± 87 MeV. Two additional 
J PC = 1— states, the V(4360) and the V(4660), have 
been reported in ISR production, but only in the reaction 
e+e" V(2S> + 7r- QUE)]. 

In this paper we present an ISR study of the reaction 
e + e~ — > J/vjit + tt~ in the center-of-mass (cm.) energy 
(E cm ) range 3.5 - 5.5 GeV. In the J/ipTT + iT~ mass re- 
gion below ~ 3.7 GeV/c 2 the signal due to the decay 
i/j(2S) — > J/ipir + ir~ dominates. A detailed comparison to 
ijj(2S) Monte Carlo (MC) simulation yields values of the 
cross section and partial width to e + e~. The high-mass 
tail of the tp(2S) MC distribution describes the data up 
to ~ 4 GeV/ c 2 quite well, and so we perform a maximum 
likelihood fit over the 3.74-5.50 GeV/c 2 mass region in 
which the fit function consists of the incoherent superpo- 
sition of a nonresonant, decreasing exponential function 
describing the J/vjtt + tt~ mass region above 3.74 GeV/c 2 
and a Breit-Wigner (BW) function describing produc- 
tion and decay of the Y(4260). Non-J/^ background is 
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41 treated by means of a simultaneous fit to the mass dis- 

42 tribution from the J/ip sideband regions. 

« This analysis uses a data sample corresponding to 

44 an integrated luminosity of 454 fb -1 , recorded by the 

45 BABAR detector at the SLAC PEP-II asymmetric-energy 

46 e+e - collider operating at cm. energies 10.58 and 10.54 

47 GeV. The detector is described in detail elsewhere [T7] , 

48 Charged-particle momenta are measured with a tracking 

49 system consisting of a five-layer, double-sided silicon ver- 
so tex tracker (SVT), and a 40-layer drift chamber (DCH), 

51 both of which are coaxial with the 1.5-T magnetic field 

52 of a superconducting solenoid. An internally reflect- 

53 ing ring-imaging Cherenkov detector, and specific ion- 

54 ization measurements from the SVT and DCH, provide 

55 charged-particle identification (PID). A CsI(Tl) electro- 

56 magnetic calorimeter (EMC) is used to detect and iden- 

57 tify photons and electrons, and muons are identified using 

58 information from the instrumented flux-return system. 

59 We reconstruct events corresponding to the reaction 

60 e + e~ — > 7isr J/ipir + ir~, where 7isr represents a photon 

61 that is radiated from the initial state e ± , thus lowering 

62 the cm. energy of the e + e~ collision which produces the 

63 J/tp /t + ?t~ system. We do not require observation of the 

64 ISR photon, since it is detectable in the EMC for only 

65 ~ 15% of the events. This is because the ISR photon is 

66 produced predominantly in a direction close to the e + e~ 

67 collision axis, and as such is most frequently outside the 

68 fiducial region of the EMC. 

69 We select events containing exactly four charged- 

70 particle tracks, and reconstruct J/ip candidates via their 

72 decay to e + e~ or For each mode, at least 

73 one of the leptons must be identified on the basis 

74 of PID information. When possible, electron candi- 

75 dates are combined with associated photons in order 

76 to recover bremsstrahlung energy loss, and so improve 

77 the J/ip momentum measurement. An e + e~ 98 

78 pair with invariant mass within (—75,-1-55) MeV/c 2 99 

79 ((—55, +55) MeV/c 2 ) of the nominal J/ip mass [T5] is ac- 100 
so cepted as a J/ip candidate. We refer to the combination 101 
ai of these e + e~ and mass intervals as "the J/ip signal 102 

82 region" . Each J/ip candidate is subjected to a geometric 103 

83 fit in which the decay vertex is constrained to the e + e~ 1M 

84 interaction region. The x 2 probability of this fit must be 105 
as greater than 0.001. An accepted J/ip candidate is kine- 106 

86 matically constrained to the nominal J/ip mass |18j and 107 

87 combined with a candidate 7r + 7r~ pair in a geometric fit 108 
as which must yield a vertex-^ 2 probability greater than 109 

89 0.001. m 

90 The value of the missing-mass-squared recoiling m 

91 against the J/ipn + 7T~ system must be in the range 112 

92 (—0.50,-1-0.75) (GeV/c 2 ) 2 in order to be consistent with 113 

93 the recoil of an ISR photon. We require also that the 114 

94 transverse component of the missing momentum be less 115 

95 than 2.25GeV/c. If the ISR photon is detected in the 116 

96 EMC, its momentum vector is added to that of the 117 

97 J/ipTr + ir~ system in calculating the missing momentum. 118 
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FIG. 1. (a) The background-subtracted data, and MC simu- 
lation modified as described in the text, for the ip(2S) peak 
region, (b) The corresponding distribution for the mass re- 
gion below 4.0 GeV/c 2 . 



The candidate t: + tt~ system has a small contamination 
due to e + e _ pairs from photon conversions. We compute 
the pair mass m c + c - with the electron mass assigned to 
each candidate pion, and remove events with m c + c - < 50 
MeV/c 2 . We estimate the remaining background by us- 
ing events that have an e + e~(/i + /i~) mass in the J/ip 
sideband (2.896,2.971) or (3.201,3.256) ((2.936,2.991) or 
(3.201, 3.256)) GeV/c 2 after satisfying the other signal re- 
gion selection criteria. 

The J/ipTt + n~ invariant-mass distribution in the re- 
gion below 4 GeV/c 2 is dominated by the ip(2S) signal. 
The peak region, after subtraction of background from 
the J/ip sideband regions, is shown in Fig. [IJa) (solid 
dots). The open dots indicate the ip(2S) MC distribu- 
tion, modified as described below. The data distribu- 
tion above ~ 3.75 GeV/c 2 (Pig.Qb)) may be due to the 
ip(2S) tail and a possible J/ip^ix^ continuum (i.e. non- 
resonant) contribution. In order to investigate this we 
performed a detailed comparison of the ip(2S) signal in 
data and in MC simulation. For the latter, we used the 
MC generator VECT0RISR [15] and a simulation of the 
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us BABAR detector based on Geant4 [30] . The resulting MC175 

120 events were subjected to the reconstruction procedures™ 

121 which were applied to the data. 177 

122 We first measured the peak mass position for both dis-™ 

123 tributions. We performed a x 2_n t of a parabola to the™ 

124 data and MC distributions in intervals of 0.5 MeV/c 2 ™ 

125 for the region within ± 5 MeV/c 2 of the nominal ip(2S)wi 

126 mass [18 . For the data, this gave a peak mass value™ 

127 of 3685.32 ± 0.02 (stat) MeV/c 2 , which is 0.77 ± 0.04™ 

128 MeV/c 2 less than the nominal value [IS]- For the MC™ 

129 events, the result was 3685.43 ± O.Ol(stat) MeV/c 2 ,™ 

130 which is 0.66 ± 0.01 MeV/c 2 smaller than the input™ 

131 value [18]. This difference is attributed to final-state- 137 

132 radiation effects. The larger deviation obtained for data™ 

133 may result from under-estimated energy-loss corrections,™ 

134 and/or magnetic field uncertainty [2"Tl |2"2"] . Each MC™ 

135 event was then displaced by 0.11 MeV/c 2 toward lower 191 

136 mass, and the parabolic fit to the new MC distribution™ 

137 was repeated. The MC distribution was normalized to™ 
us the data by using the data-to-MC ratio of the maxima™ 

139 of the fitted functions. In order to improve the MC-™ 

140 data resolution agreement, a x 2 function incorporating™ 

141 the data-MC histogram differences and their uncertain-™ 

142 ties was created for the region within ± 10 MeV/c 2 of the™ 

143 peak mass value. In the minimization procedure eacli2oo 

144 MC event was represented in mass by a superpositiomoi 

145 of two Gaussian functions with a common center, but202 

146 different fractional contributions, and normalized to one203 

147 event. The root-mean-squared (r.m.s.) deviations of the204 
us Gaussian functions, and the fractional contribution of the2os 

149 narrower Gaussian function to the normalized distribu-206 

150 tion, were allowed to vary in the fit, and the contribution207 

151 of each smeared MC event to each histogram intervahoa 

152 was accumulated. This procedure yielded a new MC his-209 

153 togram to be used in the fit to the data histogram. We2io 

154 iterated the above procedure until the change in \ 2 wasau 

155 less than 0.1, at which point the narrow (broad) Gaussian2i2 

156 r.m.s. deviation was 0.7 (6.3) MeV/c 2 , and the fractionaho 

157 contribution was 0.88 (0.12). 214 

158 In Fig. [ija) the final MC distribution is compared t02is 

159 the data in the fit region, and the agreement is good2i6 

160 (x 2 /NDF = 30.7/35, probability = 67.6 %; NDF is the 2 i7 

161 number of degrees of freedom). We integrate this MC218 

162 distribution over the entire lineshape in order to esti-219 

163 mate the VK 2 ^) signal yield in data, and obtain 20893220 

164 ± 145 (stat) events. We use the efficiency and the22i 

165 distributed luminosity (obtained from the nominal inte-222 

166 grated luminosity and the second-order radiator function223 

167 from Ref. [23 ) to obtain the cross section value 14.05 ±224 

168 0.26 (stat) pb for radiative return to the ip(2S). This is225 

169 in agreement with a previous measurement [14 . In ad-226 

170 dition we extract T(tp(2S) e+e~) = 2.31 ± 0.05 (stat)^ 

171 keV, in excellent agreement with Ref. [TSJ. 22a 

172 In Fig. [ijb) we compare the modified VK 2 ^) MC dis- 

173 tribution to the data in the region below 4.0 GeV/c 2 . The 

174 MC low-mass tail is systematically below the data distri- 



bution, but the high-mass tail provides a good descrip- 
tion of the observed events. However, we note that the 
extrapolation to this region requires the use of the ^(2S) 
Breit-Wigner lineshape at mass values which are as much 
as 1000 full-widths beyond the central mass. The exis- 
tence of many other measured final state contributions 
to the J PC = 1 amplitude in this mass region must 
call this procedure into question. Although our model 
adequately describes the data between the ip(2S) peak 
and ~ 4.0 GeV/c 2 , we cannot discount the possibility 
of a contribution from an e + e~ — > J/tjjn + Tr~ continuum 
cross section in this region. In this regard, the failure of 
the MC lineshape to describe the data in the region of 
the low-mass tail might be due to the threshold rise of 
just such a continuum cross section. 

The J/i/>7r + 7r~ mass distribution corresponding to the 
J/ip signal region is shown from 3.74 to 5.5 GeV/c 2 in 
Fig. 2(a). The shaded histogram, which has been ob- 
tained by linear interpolation from the J ftp sideband re- 
gions, represents the estimated background contribution 
to the J/ip signal region. The signal distribution shows 
an excess of events over background above 3.74 GeV/c 2 
which might result from the ip(2S) tail and a possible 
J/tjjTr + Tr~ continuum contribution, as discussed with re- 
spect to Fig. [ljb) . At higher mass we observe clear pro- 
duction of the Y(4260), and beyond ~ 4.8 GeV/c 2 the 
data are consistent with background only. There is a 
small excess of events near 4.5 GeV/c 2 , which we choose 
to attribute to statistical fluctuation. In this regard, we 
note that no corresponding excess is observed in Ref. [14] . 
The background contribution is featureless throughout 
the mass region being considered. 

In order to extract the parameter values of the 
Y(4260), we perform an unbinned, extended-maximum- 
likelihood fit in the region 3.74-5.5 GeV/c 2 to the 
J/ipir + ir~ distribution from the J/ip signal region, and 
simultaneously to the background distribution from the 
J/if) sidebands. The background is fitted using a third- 
order polynomial in J/07r + 7r~ mass, m. The mass- 
dependence of the signal function is given by f(m) — 
e(m) ■ C(m) ■ ct(to), where e(m) is the mass-dependent 
signal-selection efficiency from MC simulation with a 
J/ipir + ir~ phase space distribution, and C{m) is the 
mass-distributed luminosity [23], where we ignore the 
small corrections due to initial-state emission of addi- 
tional soft photons; e(m) increases from 9.5% at 3.74 
GeV/c 2 to 15.5% at 5.5 GeV/c 2 , and C(m) from 35 

pb" 1 / 20 McV t0 61.3 pb" 1 / 20 MeV over tne same range. 
The cross section, er(m), is given by the incoherent sum 
<r(m) = ctnr(w) + ffBw(m), where <tnr(to) is an expo- 
nential function which provides an empirical description 
of the tp(2S) tail and possible continuum contributions; 
(TBw( m ) is the cross section for the production of the 
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4.6 4.8 



TABLE I. Systematic uncertainty estimates for the Y(4260) 
parameter values. 



m(JA|«tV)(GeV/c 2 ) 



Source 


r e+e - ■ B(%) 


Mass (MeV/c 2 ) 


T ( MeV) 


Fit procedure 


+ 1.5 

-0.5 


+0 

-l 


+2 

-l 


jVIass Scale 




ztU.O 




Mass resolution 






±1.5 


MC dipion model 


±3.6 






Decay angular 








momentum 


±3.6 


±3.5 


±7 


Luminosity, etc. 


±5.4 






(see text) 
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237 
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239 
240 
241 
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243 
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248 
249 
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- 254 



FIG. 2. (a) The J/ipn + Tr~ mass spectrum from 3.74 GeV/c ; 
to 5.5 GeV/c 2 ; the points represent the data and the shaded 25 ^ 
histogram is the background from the J/ip sidebands; the solid256 
curve represents the fit result, and the dashed curve results 257 
from the simultaneous fit to the background; (b) the measured 258 
e + e~ — > J/tj)-K + -K~ cross section as a function of cm. energy; 
the solid curve results from the fit shown in (a). 



Y(4260), and is given by 



12nC 

2~ 



in- 



PS{m) 
PS {m Y ) 



B{J/tpTT+TT-) -ml 



(1), 



(m Y — m 2 ) 2 ± 



lyTy 



where my and Ty are the mass and width of the269 
Y(4260), r e+e - is the partial width for Y(4260) -> e+e~,27o 
B(J/?pir + Tr~) is the branching fraction for Y(4260) — 
J/^tt+tt", and C = 0.3894 x 10 9 GeV 2 pb. The func-272 
tion PS(m) represents the mass dependence of J/ipTr + Tt~n3 
phase space, and PS(my) is its value at the mass of the274 
Y(4260). In the likelihood function, (JBwim) is multi-275 
plied by B(J/ip — > l + l~), the branching fraction suni276 



of the 



and 



decay modes [18], since the fhW 



is to the observed events. In the fit procedure f(m) isne 
convolved with a Gaussian resolution function obtained2so 
from MC simulation. This function has a r.m.s. devia-281 



tion which increases linearly from 2.1 MeV/ c 2 at ~ 3.5 
GeV/c 2 to 5MeV/c 2 at - 4.3 GeV/c 2 . The results of the 
fit are shown in Fig.[2]Ja). The parameter values obtained 
for the Y(4260) are my = 4244 ± 5 (stat) MeV/c 2 , Ty = 
IU±H (stat) MeV, and r e+e - x B( J/^tt+tt") =9.2±0.8 
(stat) eV. 

For each J/ipTt+n^ mass interval, i, we calculate the 
e + e~ — > J/ipir + TT~ cross section after background sub- 
traction using 



^obs 



bkg 



£i ■ A • B(J/i/> -> l+l- 



(2) 



with n° bs and n^ kg the number of observed and back- 
ground events, respectively, for this interval; e.;, and Li 
are the values of e(m) and C(m) 23^ at the center of in- 
terval i. 

The resulting cross section is shown in Fig. [2jh) , where 
the solid curve is obtained from the simultaneous likeli- 
hood fit. The corresponding estimates of systematic un- 
certainty are due to luminosity (1%), tracking (5.1%), 
B(J/ip ^ 1+1-) (0.7%), efficiency (1%) and PID (1%); 
combined in quadrature. These yield a net systematic 
uncertainty of 5.4%, as indicated in Table [i]. 

The reaction e + e~ — > J/ipir + ir~ has been studied at 
the cm. energy of the ^(3770) by the CLEO [Mj and 
BES [25] collaborations. The former reported the value 
12.1 ± 2.2 pb for the e+e" ->■ ^(3770) J/^tt + tt~ 
cross section, after subtraction of the contribution result- 
ing from radiative return to the ip(2S). The dependence 
on E cm of our fitted cross section, shown by the curve in 
Fig. ||b), yields the value 31 ± 5 (stat) ± 2 (syst) pb at 
the -0(3770) with no subtraction of a ip(2S) contribution. 
This is compatible with the much more precise CLEO re- 
sult obtained after subtraction. No cross section value is 
reported in Ref . [5S] , but the results of the BES analysis 
agree within their significantly larger uncertainties with 
those from CLEO. 

The systematic uncertainties on the measured values 
of the Y(4260) parameters include contributions from the 
fitting procedure (evaluated by changing the fit range and 
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the background parametrization) , the uncertainty in the 
mass scale, the mass-resolution function, and the change 
in efficiency when the dipion distribution is simulated us- 
ing the solid histogram in Fig. 3(c), which is described 
below. In Eq. ([I]) it is assumed that Y(4260) decay to a 
J I ip and a scalar dipion occurs in an S- wave orbital angu- 
lar momentum state. However, a D-wave decay between 
the J ftp and the 7r+7r~ system can occur also, and for this 
hypothesis the fitted central values of mass, width, and 
T e+e - x B(J/ipTT+ii-) become 4237 MeV/c 2 , 100 MeV, 
and 8.5 eV, respectively. We assign half the change in 
central value of each quantity as a conservative estimate 
of systematic uncertainty associated with the decay angu- 
lar momentum. Uncertainties associated with luminos- 
ity, tracking, B(J/ip — > efficiency and PID affect 
only r e + e - • B, and their net contribution is 5.4%, as we 
discussed previously. Our estimates of systematic uncer- 
tainty are summarized in Table |TJ and are combined in 
quadrature to obtain the values which we quote for the 
Y(4260). 

We now consider the it + it~ system from Y(4260) 
decay to J/^Tr+Tr". Since the Y(4260) has I(J PC ) = 
0(1 ) and its width indicates strong decay, the tt + it~ 
system has I(J PC ) = 0(0++) or I(J PC ) = 0(2++). For 
the region 4.15 < m(J/fy> 7r+7r~) < 4.45 GeV/c 2 , the 
7r+7r~ mass distribution after subtraction of that from 
the J/i/j sideband regions is shown in Fig. [3][a). The re- 
gion below 0.32 GeV/c 2 is excluded since it is severely 
depopulated by the procedure used to remove e+e~ pair 
contamination. The distribution decreases from thresh- 
old to near zero at ~ 0.6 GeV/c 2 , rises steadily to a max- 
imum at ~ 0.95 GeV/c 2 , decreases rapidly to near zero 
again at ~ 1 GeV/c 2 , and increases thereafter. The distri- 
bution is consistent with previous measurements • 

We define 0^ as the angle between the 7r+ direction and 
that of the recoil J/tp , both in the dipion rest frame. The 
distribution in cos0,r is shown in Fig. [3]jb). The fitted 
line represents S-w&ve decay, and provides an adequate 
description of the data (x 2 / NDF = 12.3/9, probabil- 
ity = 19.7%); there is no need for a D-wave contribution, 
e.g., from / 2 (1270) — > tt + tt~ decay. 

The mass distribution near 1 GeV/c 2 suggests coher- 
ent addition of a nonresonant 7r+7r~ amplitude and a 
resonant amplitude describing the /o(980). If the peak 
near 950 MeV/c 2 is attributed to a nonresonant ampli- 
tude with phase near 90°, the coherent addition of the 
resonant / (980) amplitude, in the context of elastic uni- 
tarity, could result in the observed behavior, which is 
similar to that of the I = tt + tt~ elastic scattering cross 
section near 1 GeV (Fig. 2, p.VII.38, of Ref. [28]). How- 
ever, we have no phase information with which to support 339 
this conjecture. 340 

The distribution in Fig. ^a.) for m n7r < 0.9 GeV/c 2 341 
is qualitatively similar to that observed for the decay 342 
T(35) -> T(lS)w+iv- 27J. There, the dipion mass dis- 343 
tribution decreases from a maximum near threshold to a 344 
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FIG. 3. (a) The background-subtracted 7r+7r~ mass distribu- 
tion for the Y(4260) signal region; the dashed vertical line is 
at the nominal /o(980) mass value [18]; (b) the corresponding 
cos9^ distribution; the fitted line is for an S'-wave description; 
(c) the result of the fit using the model of Eq. Q. 



significantly non-zero minimum at ~ 0.6-0.7 GeV/c 2 , be- 
fore rising steeply toward 0.8 GeV/c 2 before being cut-off 
by the kinematic limit (0.895 GeV/c 2 ). The CLEO data 
are well-described in terms of a QCD multipole expan- 
sion [211 US] U P to m^Tr ~0.7 GeV/c 2 , but the sharp rise 
thereafter is not well-accommodated. This shortcoming 
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is more readily apparent for the much larger BABAR data4oo 
sample for this same process [30] , There the distribu-4oi 
tion begins a rapid rise toward the /o(980) region, as402 
in Fig. [3]Ja), but turns over at ~ 0.85 because of the4os 
kinematic limit at 0.895 GeV/c 2 . The CLEO multi-404 
pole expansion fit involves two amplitudes whose rela-405 
tive phase (± 155 degrees) causes destructive interfere 
ence, and hence the minimum in the mass distribution407 
at ^0.6-0.7 GeV/c 2 . The amplitudes are of similar mag-4os 
nitude in this region, and so a relative phase of =b I8O409 
degrees could yield near-zero intensity, as observed imio 
Fig. [3][a). This phase value would result in an approx-411 
imately real amplitude. However it would contain no4i2 
explicit /o(980) contribution, which seems necessary to a4i3 
description of the data of Fig.[3]Ja), and so we attempt t04i4 
describe the entire distribution using the following sim-415 
pie model. 4ie 
The nonresonant intensity distribution requires three4i7 
turning points, as in the CLEO multipole expansion de-4is 
scription, and so we choose to represent it by a fourth-419 
order polynomial, T(m 7r7r ), where the invariant 

mass of the tt + it~ system. From the phase requirement 
discussed above, it follows that the corresponding am- 
plitude can be chosen to be real and represented by 

^JTim^^). To this amplitude we add the complex S-™ 
1 421 
wave 7r + 7r~ amplitude obtained from the BABAR anal- 

422 

ysis of Dj — > ir + ir n + decay [21], which shows clear 423 
resonant behavior at the /q(980). We perform a x 2-n t424 
to the data of Fig. [3ja) using 425 



f(m n7r ) = \y / T(m^) + e i0 F /o(98O )(m 7rw )| 2 • p ■ q , (3) 427 

where .F^ggo) is proportional to the complex n + ir ~ 42g 
amplitude of Ref. |31| . and the phase <f> is determined 430 
by the fit; p is the n + momentum in the ir + n~ rest 43 i 
frame, and q is the J/tp momentum in the J/07r + 7r~432 
rest frame. We use the fitted Y(4260) mass value 433 
in calculating q, which implies a kinematic limit of 434 
1.15 GeV/c 2 for the fit function. The result is shown 435 
in Fig. Ifc). The fit is good { X 2 /NDF = 33.6/35,™ 
probability = 53.6%), and the interference contribu- 438 
tion is important for the description of the region near 439 
1 GeV/c 2 ((f) = 28° ± 24°). The / (980) amplitude44 
squared gives 0.17 ± 0.13 (stat) for the branching ratio 441 
B(J/^/o(980),/ (980) 7r + 7r")/S(J/^7r+7r-). This is 442 
somewhat smaller than the prediction of Ref. |52j , where 443 
it is proposed that the /o(98Q) contribution should be 445 
dominant. 446 
In summary, we have used ISR events to study4 47 
the process e + e~ — > J/ifjir + TT~ in the cm. en-448 
ergy range 3.74-5.50 GcV. For the Y(4260) we ob- 449 
tain my = 4244 ± 5 (stat) ±4 (syst) MeV/c 2 , 450 

T Y = 114 tit ( stat ) ±7 ( s y st ) MeV > and T e+e- x 451 
B{J/ipTT+TT-) = 9.2 ± 0.8 (stat) ±0.7 (syst) eV. These*" 
results represent an improvement in statistical precision 454 
of ~ 30% over the previous BABAR results [B], and agree 455 
very well in magnitude and statistical precision with the456 



results of the Belle fit which uses a single Breit-Wigner 
resonance to describe the data [12]. We do not con- 
firm the broad enhancement at 4.01GeV/c 2 reported in 
Ref. [H]. The dipion system for the Y(4260) decay is 
in a predominantly S'-wave state. The mass distribution 
exhibits an / (980) signal, for which a simple model indi- 
cates a branching ratio with respect to J/ipir + ir~ of 0.17 
± 0.13 (stat). 
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